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Abstract 

An analytical expression is derived for the conductance modulation of a ballistic 
two-dimensional Datta-Das Spin Field Effect Transistor (SPINFET) as a function of 
gate voltage. Using this expression, we show that the recently observed conductance 
modulation in a two-dimensional SPINFET structure does not match the theoreti- 
cally expected result very well. This calls into question the claimed demonstration 
of the SPINFET and underscores the need for further careful investigation. 
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Even two decades after the original proposal of the Datta-Das Spin Field 
Effect Transistor (SPINFET) [1], the exact analytical expression for the chan- 
nel conductance of a two-dimensional device structure has remained some- 
what obscure [see the Note at the end]. Although Datta and Das proposed 
a two-dimensional transistor structure in their original work [1], the expres- 
sion they derived for the channel conductance as a function of gate voltage 
was based on the assumption that the carrier's wavevector component trans- 
verse to the direction of current flow is zero, which effectively corresponds to 
a one-dimensional structure. No expression was derived for the conductance 
modulation in a two-dimensional structure, possibly because Datta and Das 
realized that the conductance modulation will be severely suppressed in a 
two-dimensional system. 

Recently, a report has appeared in the literature claiming demonstration of 
the Datta-Das SPINFET for the first time. The claim is predicated on the fact 
that a conductance modulation was observed in a two-dimensional SPINFET 
structure as a function of gate voltage, which could be fitted exactly with the 
equation 



where a [Vg] is the gate-controlled Rashba spin-orbit interaction strength in 
the two-dimensional channel, Vg is the gate voltage, m* is the charge carrier's 
effective mass, L is the source-to-drain separation (channel length) and is 
an arbitrary phase shift. The authors of [2] measured the expected amplitude 
A and the quantity a [Vg] in their structure independently, and then using 
as the only fitting parameter, they could fit the experimentally observed 
conductance modulation AG in their structure with Equation (1). This "fit" 
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(among others) was offered as proof that the Datta-Das transistor has been 
demonstrated. 

Ref. |2] took Equation (1) from ref. ^J, not realizing that it apphes only to a 
strictly one-dimensional channel since ref. [T] had derived it assuming that the 
wavevector component transverse to the direction of current flow is exactly 
zercQ. Equation (1) does not hold for a two-dimensional channel since there 
the transverse wavevector component will not be zero. Recently, one of us 
pointed this out [3] and derived the correct equation for a two-dimensional 
channel (of finite width) assuming that only the electron energy is conserved 
in ballistic transport. Subsequently, it was pointed out [4] that if the width 
of the channel is semi-infinite so that periodic boundary conditions can be 
imposed along the width, then the transverse wavevector (perpendicular to the 
direction of current flow) is also a good quantum number and will be conserved 
in ballistic transport. This is reminiscent of two-dimensional coherent resonant 
tunneling devices of semi-infinite width, where the transverse wavevector is 
conserved during tunneling [S]. 

Conservation of the transverse wavevector greatly simplifies the equation de- 
rived in [3]. Additionally, low temperature and low bias conditions cause fur- 
ther simplification, resulting in the following simple equation for the conduc- 
tance modulation in a two-dimensional SPINFET: 



^ The expression derived in ref. [T] did not contain the phase shift cj), but it could 
arise if a [Vq] / when Vq = 0. 




(2) 



3 



where S is a constant and 



e{kF,h,a[VG]) 



{2m* a [Vg] I^- 



(3) 




Here, is the transverse wavevector component (along the width) and /ci? is 
the Fermi wavevector. Derivation of Equation (2) is given in Appendix I. 

Clearly, Equation (2) has no similarity with Equation (1). Therefore, the con- 
ductance modulations in the one- and the two-dimensional cases are very 
different. Particularly, Equation (1) predicts that the conductance modula- 
tion could reach 100%, whereas Equation (2) shows unambiguously that it 
will never reach 100% because ensmeble averaging represented by the integra- 
tion over the transverse wavevector component k-^ will dilute the modulation 
considerably. Only in a strictly one-dimensional channel where Equation (1) 
holds, the conductance modulation can be 100%, while in a two-dimensional 
channel, it will never be 100%. 

Ref. [6] has independently derived Equation (2) for a two-dimensional SPIN- 
FET and found that it can be approximated as 



Equation (4) does not quite match Equation (1) either since the amplitude 
of the cosine function in Equation (4) is not constant, but gate- voltage de- 
pendent. Therefore, Equation (2) or Equation (4) cannot be reconciled with 
Equation (1). However, ref. [6] also found that for the particular experimental 
parameters of ref. [2], Equation (1) and Equation (2) yield similar curves for 
AG versus Vq over the range of Vq used in the experiment. This similarity 
is coincidental and will not be sustained over extended ranges of Vq- More 



AG 



cos 2m*a[VG\L/h^ + TX/A . 



(4) 
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importantly, we have found that if we use the values of m*, a [Vg], and 
L reported in [2], then the AG versus Vq curve computed from the correct 
Equation (2) does not match the experimental AG versus Vg curve reported in 
[2] very well. We show these two curves in Fig. 1. This disagreement between 
the correct theoretical result and the experimental observation casts doubt on 
the claimed demonstration of the Datta-Das SPINFET. 

We emphasize that the above disagreement however does not establish con- 
clusively that the Datta-Das SPINFET was not demonstrated in [5]. Instead, 
it casts doubt on the claimed demonstration and highlights the need for fur- 
ther investigation. Finally, the important question is if the observed voltage 
modulation was not due to the Datta-Das effect, what could it have been due 
to? Ref. [2] showed that the conductance modulation AG versus Vq disap- 
peared if the source and drain contacts were magnetized in a direction such 
that they injected and detected spins parallel to the effective magnetic field 
caused by the Rashba interaction. The modulation reappeared when the di- 
rection of magnetization was rotated by 90° so that the injected spins became 
perpendicular to the effective magnetic field. This is consistent with the Datta- 
Das effect which relies on precession of the injected spins around the effective 
magnetic field caused by Rashba interaction. Since precession cannot occur 
if the spins are parallel to the effective magnetic field, the Datta-Das mod- 
ulation will disappear in that case. The precession will occur if the injected 
spins are perpendicular to the effective magnetic field, so that the Datta-Das 
effect is recovered when the contacts' magnetizations are rotated by 90°. This 
observation is certainly supportive of the Datta-Das effect, but it could also 
be caused by other phenomena. One likely phenomenon is Ramsauer reso- 
nances in the channel [7] which can also give rise to a voltage modulation AG 
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versus Vg- Ramsauer resonances are exacerbated by a magnetic field in the 
direction of current flow. In tlie experiment, wfien tlie contacts were magne- 
tized in tlie direction perpendicular to the effective magnetic field caused by 
the Rashba interaction, they caused a real magnetic field to appear in the 
channel in the direction of current flow. This could have induced Ramsauer 
resonances. When the direction of magnetization of the contacts was rotated 
by 90°, the real magnetic field in the channel disappeared, which could have 
quenched or abated the Ramsauer resonances. Thus, the observed effect is also 
consistent with Ramsauer resonances. Consequently, further tests are required 
to identifythe origin of the observed conductance modulation unambiguously. 
The expected oscillation periods for Ramsauer resonances and the Datta-Das 
effect are of course very different, but since barely one oscillation period was 
observed in the experiment of ref. [2], it is difficult to discriminate between 
these two effects from the observed modulation. 

In summary, we have shown that the conductance modulation observed in ref. 
[2] cannot be fitted very well by the correct equation governing such a device, 
contrary to the claim of ref. P]. Moreover, there can be alternate explanations 
for the origin of the observed conductance modulation of the device. Therefore, 
further careful study is required to resolve these controversies definitively. 

Note: After the submission and acceptance of this work, we became aware of 
a paper [M. G. Pala, M. Governale, J. Konig and U. Ziilicke, Europhys. Lett., 
65, 850 (2004)] which has derived an analytical expression for the channel 
conductance of a two-dimensional SPINFET as a function of different orienta- 
tions of the contacts' magnetization. That expression reduces to Equation (2) 
when the contacts are magnetized in the +x-direction. We thank Prof. Ulrich 
Ziilicke for bringing this to our attention. 
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Appendix I 



In this Appendix, we derive the expression for the channel conductance of a 
two-dimensional Datta-Das SPINFET as a function of gate voltage. 

Consider the two-dimensional channel of a Spin Field Effect Transistor (SPIN- 
FET) in the x-z plane (shown in Fig. 2(a)), with current flowing in the 
x-direction. An electron's wavevector components in the channel are desig- 
nated as kx and k^, while the total wavevector is designated as kf. Note that 
=^1 + ^1 shown in Fig. 2(b). 

The gate terminal induces an electric field in the y-direction which causes 
Rashba interaction. The Hamiltonian operator describing an electron in the 
channel is 

2 2 

H = [I] + « [Vg] {a,px - a^p.) , (5) 



where the p-s are the momentum operators, the as are the Pauli spin matrices 
and [I] is the 2x2 identity matrix. Since this Hamiltonian is invariant in both 
X- and z-coordinates, the wavefunctions in the channel are plane wave states 
^i{ka:x+k^z) _ Consequently, in the basis of these states, the Hamiltonian is 



H 



2m 



l + a [Vg] kx -a [Vg] k^ 



-a [Vg] k^ |^ - « [Vg] k 



(6) 



Diagonalization of this Hamiltonian yields the eigenenergies and the eigen- 
spinors in the two spin-split bands in the two-dimensional channel: 
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El = - a [Vg] h (lower band); = + a [Vg] h (upper band). 

2m* 2m* 

(7) 



in 



sinO 
cos9 



flower band); 



—cos9 
sinO 



(upper band) 



where 9 = {l/2)arctan{kz/kx). The energy dispersion relations in the two 
bands (one broken and the other sohd) are plotted in Fig. 3. Note that an 
electron of energy E has two different wavevectors in the two bands given by 
k\^^ and kf\ 

We will assume that the source contact of the SPINFET is polarized in the 
+x-direction and injects +x-polarized spins into the channel under a source- 
to-drain bias. We also assume that the spin injection efficiency at the source is 
100%, so that only +x-polarized spins are injected at the complete exclusion of 
— a;-polarized spins. An injected spin will couple into the two spin eigenstates 
in the channel. It is as if the x-polarized beam splits into two beams, each 
corresponding to one of the channel eigenspinors. This will yield: 



1 

71 



+x — polarized 



1 




sinO 




—cosQ 




= Ci 








1 




cosO 




sind 



(9) 



where the coupling coefficients Ci and C2 are found by solving Equation ([9]). 
The result is 



Ci = Ci (k^, kz) = sin {0 + 7r/4) 

C2 = C2{k:,,k,) = -COS (^ + 7r/4) 



(10) 
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Note that the couphng coefficients depend on and k^. 



At the drain end, the two beams recombine and interfere to yield the spinor of 
the electron impinging on the drain. Here, we are neglecting multiple reflection 
effects between the source and drain contacts in the spirit of ref. [1]. Since 
the two beams have the same energy E and transverse wavevector k^ (these 
are good quantum numbers in ballistic transport), they must have different 
longitudinal wavevectors k^''' and k^'' since k^^'^ ^ kf'\ Therefore, these two 
beams have slightly different directions of propagation in the channel. In other 
words, the channel behaves like a "birefrigent" medium where waves with anti- 
parallel spin polarizations travel in slightly different directions. 

Hence, the spinor at the drain end will be: 



drain 



sinO 
cosO 



L+k^W 



+ C2 



-cosd 



sinO 



i \ kpL+k^W 



Jk^W 



sin {6 + n/A) sin6 



sin (6* + 7r/4) cost 



iki^''L 



COS (6* + 7r/4) sinOe 



in) 



where L is the channel length (distance between source and drain contacts) and 
W is the transverse displacement of the electron as it traverses the channel. 



Since the drain is polarized in the same orientation as the source, it transmits 
only +x-polarized spins, so that spin ffitering at the drain will yield a trans- 
mission probability |Tp where T is the projection of the impinging spinor on 
the eigenspinor of the drain. It is given by 
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V2 



1 1 



sin {6 + 7r/4) sinOe^''^^^^^ + cos {6 + 7r/4) cosOe'-^'^^^ ^ 



sin [6 + 7r/4) cosOe'-^'^^^ ^ — cos [9 + 7r/4) sinOe'-^'^^^ ^ 
sin' {6 + 7r/4) e^'^^"^ + cos^ + 7r/4) e*^^" 



Here, we have assumed 100% spin filtering efficiency. 



Therefore, 



Jk^W 



(12) 



\T\^ = cos^ {d + u/A) 



l + tan^ (^ + 7r/4)e 



: cos^ {9 + 71/4) + sm^ + 7r/4) + -cos2(2^)cos(eL), 

2 



where 6 = fcW - A;!^). 



(13) 



From Equation ([8]), we get that k^^'^ — kf''' = —2m*a \Vg\ /f?- Expressing the 



wavevectors in terms of their x- and z-components, we get: 



k'i^ ' + kl-\l klf' ' + kl = -2m* a 



,(2) 



(14) 



which yields 



-2m*a [^g] jh^kf^ + 2 (m*)' [Vg] /^^ 



X X 



, l(2) 



/2 



(15) 



From Equation (IHl), we also get that 



,(2) _ m*a [Vg] 



f9 



± 



' m*a [V( 



G\ 



+ /eg ^ fco + 



3m*a[1/G] 
2 ^2 



(16) 



where fco = \/2m*E/h. 



Now, if a [Vg] is small, then fc^^^ + k^^^ /2 ^ yko — k^. Substituting these 
results in Equation f |T5l) . we get 
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{2m*a [Vg] /h^ ) h - {m*Y [Vq] /^^ 

i.2 _ U2 

{2m*a [Vg] /h^) V2^^/h - {m*f [Vg] /h^ 
2m*E/h'^ - kl 



(17) 



The current density in the channel of the SPINFET (assuming balhstic trans- 
port) is given by the Tsu-Esaki formula: 



where q is the electronic charge, Wy is the thickness of the channel (in the 
y-direction), Vsd is the source-to-drain bias voltage and /(r^) is the electron 
occupation probability at energy t] in the contacts. Since the contacts are at 
local thermodynamic equilibrium, these probabilities are given by the Fermi- 
Dirac factor. 

In the linear response regime when Vsd 0, the above expression reduces to 



SD 










dE 



J 



This yields that the channel conductance G is 



(19) 



Vsd Vsd 



nh 







dm 

dE 



(20) 



where Isd is the source-to-drain current and Wz is the channel width. 
Using Equations f lT3|) and f lT7|l . we finally get that the channel conductance is 



G = Go+^^—^ I dE I dkzCos\2e)cos{QL) 
2Txh J J 





dE 

dm 

dE 
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2m*E 



dkz 



cos{QL) 



df{E) 
dE 



(21) 



where Go is a constant independent of and hence the gate voltage. It is easy 
to show that Go = 4^ /o°° dE J dk, fl + r_MM" 



2m* E 



dE 



If the temperature is low so that ^ S{E — Ep), then the last equation 

reduces to 



— Go 



2TTh 



dk. 





kp 



cos{QL)5{E - Ef) 



Trh 



kp 



dkz 



dk^ 



1 - 



kp 



2m*E 

cos [Q{kF,k,,a[VG])L] 



^1 
k^ 



F{a[VG],L,kF,k,), 



(22) 



where 



F {a [Vg] , L, kF, k,) = ^cos' [6 {kp, k,, a [Vg]) L] /2 + ^sm^ [6 {kp, k,, a [Vg]) L] /2 | , (23) 



and kp is the Fermi wavevector. Therefore, 

k' 



AG = G-Go 



27r/i 



dk^ 



1 - 



cos[Q{kp,k,,a[VG])L\. (24) 



The last equation is identical with Equation (2). 
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Appendix II 



In this appendix, we will derive an expression for AG assuming non-ideal 
spin injection and detection. Let us call the spin injection efficiency at the 
source contact rjs and the spin ffitering efficiency in the drain contact rjD. If 
these efficiencies are less than 100%, then the probability of a +a;-polarized 
spin being injected by the source is (1 + rjs) /2 and the probability of it being 
ffitered at the drain is {1 + rj^,) /2 when both contacts are magnetized in the 
+x-direction. Therefore the contribution to AG arising from +x-polarized 
injection and +x-polarized detection is 



[AGJ_,,_^, - - dk 





kp r ,2 



klh 



COS [0 {kF,k^,a[VG\)L\ 
(25) 



Next consider the situation when the source injects a +x-polarized spin, but 
it transmits into the -x-polarized band in the drain because spin ffitering is 
imperfect. 

In this case, since the spin injected from the source is +x-polarized. Equation 
f llip is still valid for the spinor of the electron impinging on the drain. However, 
we now have to re-calculate the projection of the impinging spinor on the -x- 
polarized state in the drain, which will give 



T — 

+x,—x I— 



1 

72 



1 -1 



sin {6 + 7r/4) sinOe''^'^^^ ^ + cos [0 + 7r/4) cosOe^^'^^^^ 
sin (6^ + 7r/4) cosOc'^'^'^^ ^ — cos {6 + tt/A) sinOe'-^'^^^ ^ 



hi'^'^ cos{2e) 



Jk^W 



(26) 
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This yields that 



IT. 



cos\2e) 



+x,—x I 



[1 - cos(eL)] = - 



A-2 



{I -COS [0 iko,h,a[VG])L]}. 

(27) 



Since the probabihty of injecting a +x-polarized spin at the source contact is 
(1 + ris) /2 and the probabihty of its transmitting into the -x-polarized band 
at the drain contact is (1 — rjr)) /2, the corresponding contribution to AG wih 
be 



[AGl 



4 2nh 



dk. 



4 



{1-cos [6 {kF,k,,a[VG])L]}. 

(28) 



Now, consider the situation when the source injects a -x-polarized spin, but it 
transmits into the +x-polarized band in the drain. 

In this case, the spin injected from the source contact is -x-polarized and we 
will have to recalculate the spinor of the electron impinging on the drain. 
Equation ([9]) will now be replaced by 



V2 



1 




sinO 




—cos 9 




= c[ 








-1 




cosO 




sind 



which yields 

C[ = sin (0-7r/4) 

C'^ = -cos {6 - 7r/4) . (30) 



Therefore, the spinor of the electron impinging on the drain is 
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[^] drain — ^1 



sinO 



cosd 



ifki'-^L+k.W. 
6 ^ ^ ~r (-^9 



— e 



ikzW 



-cos9 



sind 



sin {9 - 7r/4) sinOe''^'^^^^ + cos {9 - 7r/4) cos9e'-^^^^^ 



sin {9 - 7r/4) cos9e'-^^^^^ - cos {9 - 7r/4) sin9e'-^^^^^ 



(31) 



The projection of this spinor on the +x-polarized state in the drain gives 



T—x,+X I— 



1 

3 



1 1 



sm {6 - 7r/4) sin9e'-''^^^^^ + cos {9 -n/A) cos9e''''^^^^ 



sin {9 — tt/4) cos9e^''^^ '-^ — cos {9 — tt/4) sin9e''^'^^ '-^ 



= ^e"'^^cos{29) 



JkzW 



(32) 



Therefore, once again, 



\T^x,-xf^'-^^[l-cos{eL)] 



1 - 



^2 



(33) 



Since the probabihty of injecting a -x-polarized spin at the source is (1 — rjs) /2 
and the probabihty of its transmitting into the +x-polarized band at the drain 
is(l + 7^D)/2, the corresponding contribution to AG is 



[AG]_ 



x,+x 



:i-Vs){l + VD)q'W, 
4 2nh 



{1-cos [e {kF,k,,a[VG])L]}. 

(34) 
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Finally, consider the situation when a -x-polarized electron is injected at the 
source and transmits into the -x-polarized band of the drain contact. 

In this case, Equation (13 ip will describe the spinor of the electron impinging 
on the drain and the projection of this spinor on the -x-polarized state in the 
drain will give 



T = 



1 

V2 



1 -1 



sin {6 — 7r/4) sinOe''^'^^^ + cos {6 — n/A) cosOe^''^^^^ 



sin {6 - 7r/4) cosOe''^'^^^ ^ - cos {6 - 7r/4) sinOe"-^^^^^ 



: e'^^^ \ sin" {9 - 7r/4) e'^^^^^ + cos^ {9 - 7r/4) e'^^'^^ 



Jk^W 



(35) 



Consequently, the contribution to AG will be 



:i - Vs) (1 - Vd) q'W, 
4 2nh 



cos [Q{kF,kz,a\VG\) L] 
(36) 



Using the Principle of Superposition, the total gate voltage dependent con- 
ductance modulation will be 



AG: 



Iirn 







1 - 



k'p 



4 4 

cos [Q {kF,k^,a[VG]) L] 



■VsVd- 



27ih 



^ ' dky 



kjp 



cos [e {kF,k^,a[VG]) L] . 



(37) 



Therefore, non-ideal spin injection and filtering reduces the amplitude of any 
non-local voltage modulation by the factor rjsrju. 
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Fig. 1. Plots of channel conductance modulation AG of a two-dimensional SPIN- 
FET versus gate voltage Vg- The solid lines are theoretical results calculated from 
Equation (2) where we have used the values of m* , kp and a [Vg] reported in ref- 
erence PJ and the points are (approximated) experimental results reported in [2]. 
The amplitudes of the theoretical plots are adjusted to match the experimental 
results as closely as possible. Note that neither the periods, nor the phases of the 
experimental plots agree very well with the theoretical plots. The plots are for two 
different channel lengths of L = 1.65 /um and 1.25 ^m used in the experiments of 
ref. [2j. The experiments were carried out at low temperatures and biasesThe spin 
splitting energy in the dot specified in Fig. 2 as a function of Rashba interaction 
strength. 
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Fig. 2. (a) A two-dimensional SPINFET channel, and (b) the wavevector compo- 
nents in the plane of the channel. 
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Fig. 3. Schematic representation of the dispersion relations in the two spin spht 
bands, under the influence of the gate voltage inducing Rashba interaction in the 
channel. 
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